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It is shown that the addition of a non-linear term to the Lagrangian of the electromagnetic 
field yields a fluid with an asymptotically super-negative equation of state, causing an accelerated 
expansion of the universe. Some general properties of nonlinear electromagnetism in cosmology are 
also discussed. 
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I. INTRODUCTION 

Observations of the redshift of supernovae type la 
and the cosmic microwave background (CMB)[2| show 
that the universe is undergoing a phase of accelerated 
expansion. The various possibilities that have been ex- 
amined to account for this acceleration can be grouped 
in two classes. In the models in the first class, matter 
with unusual properties is added to the right-hand side 
of Einstein equations. Specifically, this matter (dubbed 
dark energy, although a more appropriate name may be 
smooth tension p]) must have a large negative pressure 
and, differently to the case of dark matter, it must not 
cluster with visible matter. Amongst the most popular 
candidates for dark energy, we can mention two. The 
first one is the cosmological constant A, which yields 
accelerated expansion due to the fact that p\ — —pA. 
However, it is strongly disfavoured by the (120 orders 
of magnitude) disagreement between the observed value 
and that naively predicted from quantum field theory 
0, which implies a rather unnatural fine tunning. 
The second candidate is the so-called quintessence [3, 
the simplest model of which is composed by a scalar field 
coupled to matter and under the influence of an ad- 
hoc self-interaction potential . The acceleration occurs 
in a regime in which the potential energy takes over the 
kinetic energy. Although quintessence has been success- 
ful in describing accelerated expansion, some problems 
remain: the quintessence field has not been identified, 
and a derivation of the potential from first principles is 
still lacking. 

The second class of models that try to account for ac- 
celerated expansion involve a modiflcation of the dynam- 
ics of the gravitational field in the low-curvature regime. 
This modiflcation can be achieved by considering (in 5- 
dimensional scenarios) the effects of the bulk on the dy- 
namics of gravitation on the brane or by directly 
adding to the 4-dimensional gravitational action terms 
with negative powers of the curvature scalar. An exam- 
ple of the latter type is given by the action ^3 



S 



where a is a new fundamental parameter with units of 
mass. In this article a third class of models will be intro- 
duced, which is a hybrid between the two classes already 



mentioned. We shall assume that the action for the elec- 
tromagnetic fleld is that of Maxwell with an extra term, 
namely 




(1) 



where F = F^^F^^^ . From the conceptual point of view, 
this phenomenological action has the advantage that it 
involves only the electromagnetic fleld, and does not in- 
voke entities that have not been observed (like scalar 
fields) and/or speculative ideas (like higher-dimensions 
and brane worlds). At high values of the field invariant 
F, the dynamics will be that of Maxwell plus corrections 
which are regulated by the parameter 7, while at low 
values of F it is the 1/F term that dominates ■ As 
shown in the Appendix, this modification of electromag- 
netism seems to be consistent with observation for small 
enough values of 7. Notice also that the Lagrangian in S 
is gauge invariant. Hence charge conservation is guaran- 
teed in this theory. 

The action given in Eqn. is actually only an example 
of a more general class of actions for the electromagnetic 
field with Lagrangians that can be written as 

C = Y.'^uF\ 



where the sum may involve both positive and negative 
powers of F. We shall see that in certain cases the "elec- 
tromagnetic fiuid" described by this Lagrangian can be 
thought of as composed by several noninteracting flu- 
ids, each with an equation of state (EOS) of the type 
Pi = WiPi, allowing for negative values of some of the Wi. 

Taking the electromagnetic (EM) fleld described by the 
action given in Eqn.(Q as source of Einstein's equations, 
a toy model for the evolution of the universe will be stud- 
ied. It will be shown that this model displays accelerated 
expansion caused by the dominance of the non-linear EM 
term over other forms of matter. Let us recall that un- 
der the parametrization p = wp for the dark energy, the 
universe will accelerate if w < —1/3. Values of w be- 
yond w = —1, which corresponds to the cosmological 
constant, yield a super-negative EOS which violates the 
dominant energy condition. However, observation allows 
for < —1 1 so it is worthwhile considering models 
with this type of EOS parameter. Models with w < —1 



2 



in the literature are usually constructed with a scalar 
field with negative kinetic energy term , the so-called 
phantom field. The evolution of a universe dominated 
by this type of matter is very peculiar, ending in a Big 
Rip : the energy density may grow without limit with 
time, and there is a singularity in the future. As we shall 
see, the non-linear electromagnetic model studied here 
yields ordinary radiation plus a dark energy component 
with w < —1. 



II. NONLINEAR ELECTROMAGNETISM IN 
COSMOLOGY 

The effects of a nonlinear electromagnetic theory in 
a cosmological setting have been studied in several arti- 
cles. The authors of showed that nonlinear correc- 
tions coming from an Euler-Heisenberg-like Lagrangian 
can be important in the very early universe, leading to 
the avoidance of the singularity. Some attention has been 
devoted to the dynamics of a universe governed by mat- 
ter described by Born-Infeld theory p6i |. In particular, it 
has been shown that under certain assumptions (as for 
instance the existence of a compactified space) a Born- 
Infeld field as a source yields accelerated expansion 0| . 

We shall study first some general properties of nonlin- 
ear electrodynamics in cosmology, and then pay attention 
to the specific case of the action given in Eqn.lQJ. Let us 
start with the theory defined by the Lagrangian 



(2) 



where the sum can have both positive and negative pow- 
ers of the field invariant. Due to the isotropy of the spa- 
tial sections of the Friedman-Robertson- Walker (FRW) 
model, an average procedure is needed if electromagnetic 
fields are to act as a source of gravity jlTj . Let us define 
first the volumetric spatial average of a quantity X at 
the time t by 



X 



lim — 

V-^Vo V 



-5 d X, 



(3) 



where V = J yj—g d?x and Vb is a sufficiently large time- 
dependent three-volume. In this notation, the electro- 
magnetic field can act as a source for the FRW model if 



= 0, B,: = 0, E,B,=0, 



EiEj — gij, BiBj — --^B gij. 



(4) 



(5) 



With these conditions, the energy-momentum tensor of 
the EM field associated to Maxwell's Lagrangian can be 
written as 



ip + p)v^v^ -p gf,u, 



(6) 



where 



B' 



Under the same assumptions, a general nonlinear La- 
grangian yields the stress-energy tensor given by Eqn.@ 
with 



p=-L- AE'^Lf, 



P = L+-{E'-2B')Lf, 



(7) 



where Lp = ^ . As we shall see below, the case in 
which E^ = is the only relevant case in cosmology. It 
is straightforward to prove from Eqn.lO that when there 
is only a magnetic field, the fluid can be thought of as 
composed of a collection of non-interacting fluids indexed 
by k, each of which obeys the equation of state 



Pk 



4/c 



1 Pk- 



(8) 



In particular, the fluid coming from Eqn.(Q is composed 
of ordinary radiation with = i and of another 
fluid with EOS 



Pi-i) 



It is precisely this component with negative pressure that 
may drive the acceleration of the universe. 

Let us remark that since we are assuming that Bi = 0, 
the magnetic field induces no directional effects in the 
sky, in accordance with the symmetries of the standard 
cosmological model. 

From Einstein's equations, the acceleration of the uni- 
verse is related to its matter content by 



3- = -i(p + 3p). 
a 2 



(9) 



In order to have an accelerated universe, the matter must 
satisfy the constraint (p -I- 3p) < 0. Assuming that the 
non-linear electromagnetic field is the dominant source 
of gravity and using the quantities defined in Equ.© we 
can write 

P + 3p = 2{L-4B'^Lf). 
Hence the constraint (p + 3p) < translates into 



Lp > 



(10) 



It follows that any nonlinear electromagnetic theory that 
satisfies this inequality yields accelerated expansion un- 
der the conditions given in Eqns.(4) and (5). 
As an aside, note that from Eqn.Q, 



P 



B')Lf. 



(11) 



Consequently, nonlinear electromagnetic theories that 
have a state for which Lp = may generate an effec- 
tive cosmological constant. 
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A. A particular example 

We shall consider next the action given by Eqn.(^ in 
the cosmological setting given by the FRW metric. The 
corresponding equations of motion (which reduce to or- 
dinary electromagnetism when 7 = 0) are 



1 + 



= 0. 



(12) 



From now on we shall assume that = 0, while there is 
a residual magnetic field characterized by Bi = 0, BiBj = 
— ^B'^gij. This assumption is consistent with the fact 
that the electric field will be screened by the charged 
primordial plasma, while the magnetic field lines will be 
frozen [13 • In this case, 



B^ 



P = 



1^1. 

Yb^' 



where we have set 7 — — /i^ llj. From the conservation 
law 

p + 3{p + p)-=0, 
a 

and setting ao = 1 we get that B = Bq/o^. Hence the 
evolution of the density with the scale factor is 



B^ 1 



j£_ 

2Bl 



(13) 



For small a it is the ordinary radiation term that domi- 
nates. The 1/F term takes over only after a = y^SoZ/i, 
and grows without bound afterwards. In fact, the curva- 
ture scalar is 



R 



4/ 



showing that there is a curvature singularity in the future 
of universes for which a ^ 00 (see Sect III B)) . 

Using the expression for the density in Eqn.Q gives 



EOM (12) the universe enters a phase of accelerated ex- 
pansion. The non-linear theory depends upon the yet 
undetermined parameter p. It is possible to get an order- 
of- magnitude estimation for this parameter as follows. 
Recent observations 19] suggest that the approximate 
current values for fi^e (the ratio of the dark energy den- 
sity to the critical density pc) and the Hubble parameter 
are Vide ~ 0.70 and h w 0.70. Assuming that the dark en- 
ergy can be described by the non-linear term in Eqn.(l), 
we get that P(-i) ~ O.TOpc- In other words. 



1.40pc 



(14) 



From this equation, we can get an upper bound for hp'^ 
if we take Bq as that associated to the CMB radiation. 
This can be obtained from — "^AT x 10^^ h^^. It 
follows that Bq 



in Ean. l(n|l we get 



^rad 

2.47 X 10~^ h^^pc- Inserting this value 



hp"^ w 3.74 X IQ-^^ gr/cm3 



(15) 



This value of hp'^ together with the results given in the 
appendix show that, in the case of a point charge, the 
corrections to Maxwell's electromagnetism coming from 
the non- linear term in Eqn.(Q are negligible. 



B. Dynamics 

The dynamics of the universe with matter density 
given by Eqn. (|13|l can be obtained qualitatively from the 
analysis of Einstein's equations We shall be interested 
only in the early regime (where radiation with p — 1 /3p 
drives the expansion), and in the late evolution (where 
the 1/F term is important). Consequently, dust-like mat- 
ter will be disregarded. Equation 



6 a4 2 52 



implies that d has a zero at a — Uc, while equation 



B'o 1 



3p^ 
2Bl 



0. 



To get a regime of accelerated expansion, we must have 



Bl 



„8 

3^«'<0, 
^0 



which implies that the universe will accelerate for a > Ac, 
with 



Bl 
3p» 



1/8 



We have shown that when the matter content of the uni- 
verse is dominated by a non-linear EM field with the 



where k = 8ttG, and e = 0, ±1 depending on the 3-section 
of the geometry, can be written as 



where 



d^ -I- V{a) 



k f B, 



(16) 



6 



Bq 



Ean. (|16|l can be interpreted as describing the one- 
dimensional motion of a particle with energy — e under 
the influence of the potential V{a). The potential is al- 
ways negative, and has a maximum at a = Uc- 
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FIG. 1: Qualitative plot of the potential for the case Vc > — 1. 
The region between ai and a2 is classically forbidden. 



The analysis of the flat and open cases shows that this 
system allows for universes that starting with a big bang 
at a = 0, expand with negative acceleration up to a = Oc, 
and then run away to infinity with positive acceleration, 
the difference between the two being the value of the 
Hubble parameter at a — a^. Note in particular that 
for e = 0, the dynamics is invariant under the change 
a — !■ 1/a. 

The closed case is more interesting due to the fact that 



27/4 



(17) 



can be ^ — 1 [20j. We shall examine three separate cases. 

Case (1): Vc < —1. This case displays the same be- 
haviour as that of the e = 1,-1 cases: the universe may 
start from a big bang, then a deaccelerated expansion 
follows up to a = Oc, and afterwards accelerated eternal 
expansion takes place. 

Case (2): Vc = —1- The same as the previous case, but 
this time H{ac) = 0, so that a = Oc is an unstable equi- 
librium point. 

Case (3): Vc > —1. Two different behaviours are allowed 
(see Fig.(^). The universe may start from a = 0, expand 
up to oi and then re-collapse. Alternatively, it may start 
from infinity, collapse to a = 02, and then re-expand to 
infinite. 



III. DISCUSSION 

The accelerated expansion suggests that the EOS pa- 
rameter w of the dark energy may be smaller that -1, 



the limit usually imposed by theoretical arguments such 
as the dominant energy condition. The confirmation of 
these findings will have a profound impact for fundamen- 
tal physics. In this article we have given a particular re- 
alization of a model that displays a time-dependent EOS, 
with a super-negative EOS parameter as a limit. Instead 
of using a negative kinetic energy term for a scalar field, 
the price to pay was to modify the action for the electro- 
magnetic field by the addition of a non-linear term that 
respects gauge invariance. As shown in the appendix, 
this modification is harmless in the case of a point charge 
[23! . The result of the modification in FRW cosmology is 
to accelerate the universe, ending in a Big Rip, indepen- 
dently of the spatial section of the geometry (except for 
one particular case in the closed universe, see Fig.(l)). 

Attempts are currently under way to reconstruct the 
EOS from observation [5j. This reconstruction implies 
the measurement of the jerk (third derivative of the scale 
factor) 21]. Although the toy model we analyzed here 
yields an EOS with a value (^(-i) — —7/3) that seems 
to be rather disfavored by observation [lj| , it would 
be interesting to see if the addition of other non-linear 
terms gives the chance of obtaining more realistic values 
for the EOS parameter. 



Appendix 

The least we should require of any modification to 
Maxwell's equations is consistence with the case of a 
static electric field generated by a point charged particle. 
The EOM for a general non- linear Lagrangian L — L{F) 
are 



(18) 



where J'^ is an external current. It follows immediately 
that d^J^ = 0, showing again that charge is conserved in 
a theory with L = L{F). In the case of the point charge, 
Ean. (|18|l reduces to 

r^Lp E{r) = const.. 

In the case of the lagrangian given in Eqn. (^) we get 



£;(r)4 ^£;(r)3 - / = 0, 



(19) 



where Q = 4 x const., and we have set ft = 1. Although 
this cubic equation could be solved exactly, we shall only 
need its solution for large and small values of r. For 
>> Q/ ji^ we get 
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while in the limit « Qj ^ 





(20) 



(21) 
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By taking derivatives of Eqn. (|19l) it can be shown that the 
function E{r) has no extrema |i2J|. Hence, the modulus 
of the electric field monotonically decreases with increas- 
ing r, from an infinite value at the origin to a constant 
(nonzero but small) value at infinity. This situation is 
akin to that in the theory defined by the action 




It was shown in lOj that the static and spherically 
symmetric solution of this theory does not approach 
Minkowski asymptotically; it tends instead to (anti)-de 



Sitter spacetime. Analogously to the gravitational field 
in [I0|, by choosing the parameter /z small enough the 
electric field of a point particle in the non-linear theory 
given by Eqn. ^ will pick up only extremely tiny correc- 
tions, given by Eans. lf^ and if^ . 
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